      SUBROUTINE ATMO62(H,RHO,ASOUND,TEMP,PRES)
C
C         'ATMO62' IS THE ATMOSPHERIC MODEL BASED ON THE U.S.
C        STANDARD ATMOSPHERE OF 1962
C
      IMPLICIT REAL*8 (A-H,O-Z)
      DIMENSION  TB(26),HB(26),BM(26),R(26),PB(26),TWTM(26)
      DATA  TB /
     1 320.65D0,  288.15D0,  216.65D0,  216.65D0,  228.65D0,  270.65D0,
     X 270.65D0,
     2 252.65D0,  180.65D0,  180.65D0,  210.65D0,  235.65D0,  260.65D0,
     X 360.65D0,
     3 560.65D0,  960.65D0,  1110.65D0,  1210.65D0,  1350.65D0,
     X 1550.65D0,
     4 1830.65D0,  2160.65D0,  2290.65D0,  2420.65D0,  2590.65D0,
     X 2700.65D0/
      DATA  HB /
     1 -5000.0D0,  0.0D0,  11000.0D0,  20000.0D0,  32000.0D0,
     X 47000.0D0,  52000.0D0,
     2 61000.0D0,  79000.0D0,  88743.0D0,  98451.0D0,  103294.0D0,
     X 108129.0D0,
     3 117776.0D0,  127394.0D0,  146541.0D0,  156071.0D0,  165571.0D0,
     X 184485.0D0,
     4 221967.0D0,  286476.0D0,  376312.0D0,  420240.0D0,  463526.0D0,
     X 548230.0D0,  630530.0D0/
      DATA  BM /
     1 -0.0065D0,  -0.0065D0,  0.0D0,  0.001D0,  0.0028D0,  0.0D0,
     X -0.0020D0,
     2 -0.0040D0,  0.0D0,  0.0030902D0,  0.005162090D0,  0.005170631D0,
     3 0.01036591D0,  0.02079434D0,  0.020891D0,  0.015739D0,
     X 0.0105263D0,
     4 0.00740193D0,  0.00533590D0,  0.00434048D0,  0.00367336D0,
     X 0.002959388D0,
     5 0.00300328D0,  0.00200699D0,  0.00133657D0,  0.0D0/
      DATA  R /
     1 0.193049997D+1, 0.122500263D+1, 0.363918856D-0, 0.880350657D-1,
     2 0.132250458D-1, 0.142753743D-2, 0.759434283D-3, 0.251091339D-3,
     3 0.200114310D-4, 0.317015171D-5, 0.497409776D-6, 0.211674012D-6,
     4 0.982986927D-7, 0.243629359D-7, 0.759159245D-8, 0.183655895D-8,
     5 0.115938455D-8, 0.804025792D-9, 0.434910767D-9, 0.156481698D-9,
     6 .358881794D-10, .650950819D-11, .312789872D-11, .157970712D-11,
     7 0.464882739D-12, 0.154055260D-12/
      DATA  PB /
     1 0.120676249D+4, 0.472681269D+3, 0.105578616D+3, 0.255403660D+2,
     2 0.404931098D+1, 0.517378375D+0, 0.275239632D-0, 0.849500768D-1,
     3 0.484093063D-2, 0.766885906D-3, 0.140309744D-3, 0.667954773D-4,
     4 0.343097486D-4, 0.117659643D-4, 0.569949880D-5, 0.236255541D-5,
     5 0.172430381D-5, 0.130345970D-5, 0.786596574D-6, 0.324927509D-6,
     6 0.879763417D-7, 0.188339643D-7, 0.959446043D-8, 0.512056269D-8,
     7 0.161272861D-8, 0.557126634D-9/
      DATA  TWTM /
     1 28.9644D0,  28.9644D0,  28.9644D0,  28.9644D0,  28.9644D0,
     X 28.9644D0,
     2 28.9644D0,  28.9644D0,  28.9644D0,  28.9644D0,  28.88D0,
     X 28.75D0,  28.56D0,
     3 28.07D0,  27.58D0,  26.92D0,  26.66D0,  26.40D0,  25.85D0,
     X 24.70D0,  22.66D0,
     4 19.94D0,  18.82D0,  17.94D0,  16.84D0,  16.17D0/
      DATA AG /9.80665D0/, AM0 /28.9644D0/, ARR /8.31432D+3/,
     .     AR /6356766.0D0/
      DATA CON1 /0.3048D0/, CON2 /1.8D0/, CON3 /0.0021156D0/,
     .     CON4 /263.83052D0/
C
C
      IF((H-2275000.0D0).LE.0.0D0) GO TO 2
      RHO=3.3D-15
      ASOUND=3413.42933D0
      TEMP=2713.95661D0
      PRES=0.126068033D-11
      GO TO 105
    2 ZZ=H*CON1
      AH=AR*ZZ/(AR+ZZ)
      DO 7 I=1,26
      IF(AH-HB(I)) 7,8,9
    9 IF(AH-HB(I+1)) 8,7,7
    8 K=I
      GO TO 10
    7 CONTINUE
   10 WS4=(AH-HB(K))
      TM=TB(K)+(BM(K)*WS4)
      CALL ARTLU1(1,AH,HB,WTM,TWTM)
      TEMP=TM*CON2*(WTM/AM0)
   12 IF(BM(K)-0.0D0) 15,16,15
   15 WS5=(TB(K)/TM)
      WS6=AG*AM0/(ARR*BM(K))
      RHO=(R(K)*WS5**(1.0D0+WS6))/R(2)
      PRES=(WS5**WS6)*CON3*PB(K)
      GO TO 17
   16 WS7=DEXP(-(AG*AM0*WS4)/(ARR*TB(K)))
      RHO=R(K)*WS7/R(2)
      PRES=PB(K)*WS7*CON3
   17 CONTINUE
      ASOUND=CON4*DSQRT(TEMP/WTM)
  105 RETURN
      END
